Hydrogen functions

Start with some imports fro Symbolic Python library:

Define some variables, radial, polar, azimuthal, time, and two frequencies:

Look at a few of the radial equations and the spherical harmonics

Notice that instead of Ylm the name is Ynm... the arguments to the function are still the quantum numbers l and m

Write the equation for the $|nlm\rangle = |100\rangle$ state. Use the sympy method .expand(func=True) to convert to the actual expression. To create this state, we combine the Radial function and the Ylm function. Make sure to set n, l, and m to the correct values. The fourth argument to R_nl is Z which we set to 1 since we are talking about a 1-proton nucleus.

The combination of R_nl and Ynm should look like the following (replace N, L, and M with the appropriate values):

R_nl(N, L, r, 1)*Ynm(L, M, theta, phi).expand(func=True)

Integrating over all space

Remember spherical coordinate integrals of function $f(r,\theta,\phi)$ over all space look like: $$\int_0^\infty\int_0^\pi\int_0^{2\pi}r^2\sin(\theta)drd\theta d\phi \,\,f(r,\theta,\phi)$$ so you alwasy need to add a factor of r**2*sin(theta) and then integrate r from 0 to infinity, theta from $0-\pi$ and phi from $0-2\pi$. As a check, you should integrate the square of the psi100 wavefunction over all space to see that it equals 1 (i.e. it is normalized)

Now do the $|210\rangle$ state:

Note, if you compare these to listed solutions (for example at http://hyperphysics.phy-astr.gsu.edu/hbase/quantum/hydwf.html#c3) you see that there are not any factors of $a_0$. This is because the R_nl function is defined in units of $a_0$. $a_0$ is the Bohr Radius: http://en.wikipedia.org/wiki/Bohr_radius

Now calculate $\langle z \rangle$:

To calculate $\langle z \rangle$ we need to convert to spherical coordinates: $z = r\cos\theta$. The terms in the following integral are the $r^2\sin\theta$ then $z$ (in spherical coords) then the wave function squared.

No surprise, the average z position of the electron in the hydrogen atom is 0.

Now for problem 13.21

find $\langle z \rangle(t)$. Use the same integral, but add a time-dependent piece to each term in the wavefunction, add them together and multiply by the complex conjugate.

We need to interpret this result. First you should show that this expression is simply a constant amplitude factor times $\cos((w2-w1)t)$, in other words $\langle z \rangle$ oscillates at frequency w2-w1.

Your assignment:

Explore other combinations of states and draw conclusions about the z behavior from the results. You may not be able to get these expressions to simplify, but the important thing is to look for the time dependence and simplify that part.

• Does $\langle z \rangle$ oscillate for any combination of two Hydrogen states $|nlm\rangle$?
• Are there restrictions on what n values give oscillating $\langle z \rangle$ expressions? (hint, to keep it simple, always let one state be the n=1 state)
• How does $\langle z \rangle$ change with different l and m values are used in the state?

Hints for interpreting your results:

• What are the relavant frequencies in your expression for $\langle z \rangle$ and why?
• Simplify one of your $\langle z \rangle$ expressions and write the time dependence in terms of the frequencies w2 and w1.